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Abstract. Selection of thermally suitable microhabitats plays an important role in the thermoregulation of ectotherms. We 
studied microhabitat preferences in two populations of the Balearic lizard, Podarcis lilfordi, that live on two coastal islets of 
Menorca (Balearic Islands, Spain: Aire and Colom). To assess the roles of different microclimatic traits driving lizards’ de¬ 
cisions, we employed a case/control design. We measured microclimate variables in locations used by lizards (cases) and in 
paired random unused locations (controls). By comparing paired cases vs controls, we prevent spatial and temporal varia¬ 
tion of environmental conditions. In order to test if microhabitat selection depends on thermoregulation requirements, we 
studied lizards in three situations: (i) cold’, when body temperatures (Tj^) of lizards were below their preferred temperature 
range (PTR), (2) ‘optimal’, when were within the PTR, and (3) ‘warm’, when exceeded the PTR. Substrate temperature 
was the main abiotic trait that determined microhabitat preferences of P. lilfordi, sometimes in conjunction with air tem¬ 
perature, while wind speed and humidity were not significant. ‘Cold’ lizards selected warmer microhabitats than the mean, 
but only on Colom islet. Notably, ‘optimal’ lizards also preferred microhabitats that were warmer than their surroundings, 
and ‘warm’ lizards did not select cold microhabitats in order to cool themselves, but rather acted randomly regarding tem¬ 
perature, wind and humidity. These results for ‘optimal’ and ‘warm’ lizards were consistent for both islets. We also studied 
the beginning and end of lizards’ diel activity. We found that lizards on Aire were active for approximately one hour more 
per day than lizards on Colom, which might be related to differences in melanism. However, activity started at an air tem¬ 
perature of approximately 17-19°C and finished at approximately the same temperature on both islets. 

Key words. Balearic lizard, behavioural thermoregulation, habitat choice, islands, Lacertidae, Menorca, microclimate, tem¬ 
perature, wind speed. 


Introduction 

Temperature affects all aspects of life and evolution of ecto¬ 
therms (PoRTNER 2002, Angilletta 2009). One impor¬ 
tant dimension of their thermal biology is thermoregula¬ 
tion (e.g., Hertz et al. 1993), that is, the ability to regulate 
body temperature (reviewed by Angilletta 2009). Lizards 
mainly thermoregulate by adjusting their activity periods 
(Hertz 1992, Adolph & Porter 1993), shuttling between 
different thermal microhabitats (Heath 1970, Bauwens 
et al. 1996), or adjusting their body posture (Bauwens et 
al. 1996). The combination of these strategies depends on 
the balance between costs and benefits (Huey & Slatkin 
1976, Blouin-Demers & Nadeau 2005). Other abiotic fac¬ 
tors, such as wind, can also have a significant effect on ther¬ 
moregulation (ScHEERS & van Damme 2002, Maia-Car- 
NEiRO et al. 2012, Logan et al. 2015, Ortega et al. 2017). 

Lacertid lizards are heliothermic thermoregulators (Ar¬ 
nold 1987, Castilla et al. 1999), and the selection of suit¬ 
able microhabitats is a key point of their thermoregulation 
(Bauwens et al. 1996, Castilla et al. 1999). Microhabi¬ 


tat selection is conditioned by a balance between the costs 
and the benefits of their thermal characteristics, the oppor¬ 
tunities for reproduction and feeding, and the avoidance 
of predators (e.g., Huey 1991, MartIn & Salvador 1997, 
Downes & Shine 1998, Downes 2001). The study of the 
abiotic factors that influence habitat selection at a micro¬ 
scale helps to understand the mechanisms of thermoregu¬ 
lation (ScHEERS & van Damme 2002, Row & Blouin-De¬ 
mers 2006). In addition, it can be important to compre¬ 
hend the response of reptiles to thermal changes in their 
habitats, since the selection of microhabitat involves physi¬ 
ological consequences that can be crucial to their demo¬ 
graphy and ecology (Huey 1991, Huey et al. 2003). 

Here we studied the influence of main abiotic variables, 
i.e., temperature, wind, and humidity, on the microhabitat 
choice of an insular lacertid lizard. Our main goal was to 
test how lizards selected microhabitats that allowed them 
to obtain body temperatures close to their optimal range. 
Secondly, we aimed to assess the relative importance of 
different abiotic traits on their microhabitat selection. We 
used the Balearic lizard, Podarcis lilfordi, as a model or- 
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ganism, comparing two populations that live on two islets 
that are situated close to each other, and differ in aspects of 
their thermal biology (Ortega et al. 2or4). 

Materials and methods 

Study species and area 

The Balearic lizard, Podarcis lilfordi (Gunther, r874), is 
an endemic lacertid of Mallorca and Menorca (Balearic Is¬ 
lands, Spain). We studied two subspecies, P. 1 . lilfordi, from 
Aire islet, and P. 1 . brauni (Muller, r927), from Colom is¬ 
let. Both islets are close to the coast of Menorca. 

Aire has a surface area of 35 ha, mostly covered with ha- 
lophytic shrubs. The average annual temperature is i6.7°C 
and average annual precipitation 625 mm (Jansa i Clar 
1979)- Food resources for lizards are limited, and preda¬ 
tion pressure is very sporadic. P. 1 . lilfordi is a small lacer¬ 
tid lizard (mean snout to vent length - SVL - for males = 
68.98 mm; mean weight = 9.75 g, mean SVL for females = 
6r.73 mm; mean weight = 6.34 g) in a melanistic popula¬ 
tion (Perez-Mellado & Salvador r988). Podarcis 1 . lil¬ 
fordi achieves a mean effectiveness of thermoregulation of 
0.81 on Aire in spring (Ortega et al. 2014). 

Colom has a surface area of 59 ha, and mean annual tem¬ 
perature and rainfall are similar to those on Aire. Its vegeta¬ 
tion is more diverse than on the latter, with abundant and 
larger shrubs and even some arboreal patches. Food resourc¬ 
es for lizards are richer on Colom than on Aire (unpubl. 
data). Podarcis 1 . brauni has a brownish-greenish coloration, 
and a significantly smaller body size than P. 1 . lilfordi (mean 
SVL for males = 67.96 mm; mean weight = 8.50 g, mean SVL 
for females = 58.68 mm; mean weight = 6.54 g). Podarcis 1 . 
brauni achieves a mean effectiveness of thermoregulation of 
0.86 on Colom in spring (Ortega et al. 2or4). 

Daily emergence and retreat 

For both islets, we recorded the occurrence of four events: 
(r) Emergence time: when we observed the first lizard leav¬ 
ing its shelter for its first basking session in the morning; 
(2) starting activity: when one of the first lizards basking 
started its activity (foraging, mating, fighting, etc.); (3) end¬ 
ing activity: when lizards stopped their daily activity and 
started their last basking session before retreating; (4) re¬ 
treat time: when we observed the last lizard retreating to its 
shelter after its final basking session (we waited 5 minutes 
in order to confirm it had indeed retreated). At these four 
events, we measured wind speed (in ms ')> relative humidity 
(RH, in %), and air temperature (TJ with a Testo* 925 digital 
thermometer, and a Kestrel® 3000 anemometer, respectively. 

Preferred temperature range (PTR) 

We studied the PTR of P lilfordi in spring in order to design 
our study of microhabitat selection. To obtain the PTR, we 


measured selected body temperatures (T ^j) of P. lilfordi on 
r2 and r3 June 2008 in a thermal gradient laboratory ex¬ 
periment. Lizards captured on Aire were transported to the 
laboratory in Sant Lluis (Menorca) and housed in individu¬ 
al opaque terraria. There, lizards were fed with mealworms 
and crickets, and water was available to them ad libitum. 
The thermal gradient was built in a glass terrarium (roo x 
60 X 60 cm), using a r50 W infrared lamp mounted above 
one of the ends, which created a gradient from 20 to 6o°C. 
Prior to the experiment, we allowed lizards one day of ac¬ 
climation to captivity, exposing them to the local natural 
conditions of light and temperature, which were similar to 
those in their habitat. Over the next two days, we measured 
a selected temperature of an individual lizard each hour 
from 08:00 to r7:oo h GMT with a Testo® 925 digital ther¬ 
mometer. We tested 20 adult lizards (ro males and females 
each), each of which was measured six times, obtaining a 
total of r20 (20 X 6) values of T^p of which we considered 
the 50% as the PTR. We released lizards in the place of their 
capture on Aire islet immediately after the experiment. 

Microhabitat selection 

We studied microhabitat selection on Aire and Colom in 
April and May of 2007, using a case/control design. We con¬ 
ducted this study during spring in order to exclude seasonal 
variations in thermoregulation, which is known to differ be¬ 
tween the two islets (Ortega et al. 2014). We noose-captured 
110 lizards on Aire and r35 on Colom during throughout their 
activity periods from 06:00 to r8:oo h GMT. We measured 
their body temperatures (Tj^) immediately after capture with a 
digital thermometer, keeping the probe in the shade. 

As we were interested in abiotic constraints and thermal 
influences of microhabitat choice, we sorted the lizards of 
both islets in three categories: (r) cold’ lizards (T^^ < PTR) that 
would need to warm up, (2) bptimaf lizards (T,^ within the 
PTR), and (3) ‘warm’ lizards (Tj^ > PTR) that would need to 
cool down. 

For each lizard, we described the microhabitat at its point 
of capture (case) as to its type of substrate (rock, soil, grass 
or bush), degree of insolation (full sun, filtered sun, or full 
shade), wind speed (ms ')> relative humidity (% RH), and sub¬ 
strate (T) and air (TJ temperatures. In order to assess abi¬ 
otic constraints in lizards’ microhabitat preferences, we also 
recorded aU these variables at four selected points (controls), 
which represent the availability of climatic conditions for a 
particular lizard at a certain time and place. We selected these 
four points at r m distance from the respective capture point in 
each compass direction (North, South, East and West). Mean 
values of these four points around each capture point were the 
control values for our microhabitat selection analysis. 

Data analysis 

Mean values are accompanied by standard errors (SE). 
Parametric statistics were performed when data followed 
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Table 1. Paired comparisons of the use of each microhabitat and 
insolation situation by lizards relative to the availability of alter¬ 
natives. D.f = 1 for aU comparisons (for details, see text). 




G 

P 

Colom Substrate 

Rock-soil 

14.826 

< 0.0001 


Rock-grass 

47.907 

< 0.0001 


Rock-bushes 

62.880 

< 0.0001 


Soil-grass 

2.823 

0.093 


Soil-bushes 

22.594 

< 0.0001 


Grass-bushes 

15.629 

< 0.0001 

Insolation 

Full sun-filtered sun 

16.304 

< 0.0001 

situation 

Full sun-shade 

18.227 

< 0.0001 


Filtered sun-shade 

0.644 

0.422 

Aire Substrate 

Rock-soil 

10.914 

0.001 


Rock-grass 

19.928 

< 0.0001 


Rock-bushes 

20.176 

< 0.0001 


Soil-grass 

5.163 

0.023 


Soil-bushes 

5.900 

0.015 


Grass-bushes 

0.057 

0.811 

Insolation 

Full sun-filtered sun 

28.799 

< 0.0001 

situation 

Full sun-shade 

49.981 

< 0.0001 


Filtered sun-shade 

3.716 

0.054 


the assumptions of normality and homogeneity of varianc¬ 
es. If these assumptions were not fulfilled, even after log- 
transformation, non-parametric equivalents were calculat¬ 
ed (SoKAL & Rohlf 1995, Crawley 2012). Analyses were 
conducted in R, version 3.1.3 (R Core Team 2015). 

The PTR was only used in order to establish the three 
categories of lizard body temperature (‘cold’, ‘optimal’ and 
‘warm’). First, we studied univariate microhabitat selection 
in the three thermoregulation situations (‘cold’, ‘optimal’, 
and ‘warm’) for both islets in order to explore the patterns 
of influence of the measured abiotic variables. Then, to add 
robustness to the repeated measuring design, we assessed 
microhabitat selection with a case/control procedure, us¬ 
ing capture points of lizards as cases and the mean values of 
the four points surrounding each capture point as controls 
(Compton et al. 2002). This case/control design is analo¬ 
gous to the use/no use of Manly et al. (2002). In order to 


analyse the paired data, we applied matched-pairs logistic 
regression or paired logistic regression (Compton et al. 
2002, Hosmer & Lemeshow 2004). We selected the model 
that best fitted to each set of data, using the multi-model 
selection function of AIC^ (Burnham & Anderson 2002) 
with the package MuMIn (Barton 2012). Coefficients 
from paired logistic regression models are interpreted in 
terms of odds ratios and not in absolute values (Keating 
& Cherry 2004). 

Results 

Body temperatures and climatic conditions 

In the Aire lizards, Tj^ did not differ between sexes (one¬ 
way ANOVA, = 0.131, p = 0.718) and age groups (one¬ 
way ANOVA, F^ = 0.091, p = 0.764), for which reason 
we pooled data for both sexes in subsequent analyses. In 
the Colom lizards, T^^ were also similar between sexes (one¬ 
way ANOVA, F^ = 0.546, p = 0.462) and age groups (one¬ 
way ANOVA, F^ = 0.019, p = 0.890). Body temperatures 
were similar for both islets (Colom: 32.84 ± o.36°C; Aire: 
32.85 ± o.24°C; Mann-Whitney U test, U = 6814, p = 0.268). 
Air temperatures at the points of capture were also similar 
on both islets (Colom: 26.29 ± o.33°C; Aire: 26.40 ± o.27°C; 
Mann-Whitney U test, U = 7753, p = 0.860), and substrate 
temperatures were marginally higher on Aire (Colom: 
28.72 ± o. 38°C; Aire: 29.71 ± o.39°C; one-way ANOVA, 

250 “ 3'^2.i, p = 0.058). Wind speed was significantly high¬ 
er on Aire (Colom: 0.29 ± 0.05 msAire: 0.61 ± 0.07 ms 
Mann-Whitney U test, U = 6042.5, p = 0.001), and rela¬ 
tive humidity was also significantly higher on Aire (Colom: 
61.71 ± 1.01%; Aire: 69.28 ± 0.92%; one-way ANOVA, F^ 

= 30.178, p < 0.0001, Table 1). 

Regarding thermal availability, that is, the mean values 
of temperatures at surrounding points, available substrate 
temperatures were higher on Aire than on Colom (Colom: 
26.24 ± o.4i°C; Aire: 27.66 ± o.39‘’C; one-way ANOVA, 
F^ = 6.210, p = 0.013). Available air temperatures were 
also significantly higher on Aire than on Colom (Colom: 
23.91 ± o.29‘’C; Aire: 25.00 ± o.24‘’C; Mann-Whitney U 
test, U = 4218.00, p = 0.004) and this was also the case for 
wind speed (Colom: 0.31 ± 0.05 msAire: 0.53 ± 0.06 ms 
Mann-Whitney U test, U = 4424.50, p = 0.006) and rela- 


Table 2. Time (GMT), wind speed (speed), relative humidity (RH), and air temperature (T^) of the four studied periods of the daily 
activity of Podarcis lilfordi on Aire and Colom islets (mean ± SE). 



Islet 

Time 

Speed (ms ') 

RH (%) 

T (°C) 

Emerging 

Colom (n=4) 

06:27±0:12 

0.70±0.32 

83.25±4.03 

17.05±0.51 


Aire (n=7) 

05:34±0:23 

1.65±0.74 

81.14±3.36 

18.06±0.59 

Starting activity 

Colom (n=4) 

06:57±0:07 

1.24±0.64 

79.75±4.37 

19.52±0.69 


Aire (n=7) 

06:11 ±0:16 

1.45±0.74 

84.00±4.81 

18.96±0.38 

Ending activity 

Colom (n=2) 

17:17±0:19 

1.02±0.15 

67.50±0.50 

20.85±0.05 


Aire (n=7) 

17:54±0:04 

2.29±0.48 

80.14±2.96 

20.73±0.58 

Retreating 

Colom (n=2) 

17:56±0:02 

0.67±0.05 

73.00±7.00 

18.00±0.40 


Aire (n=7) 

18:42±0:06 

2.21±0.56 

90.71±3.85 

18.51±0.35 
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tive humidity (Colom: 64.53 ± 0.89%; Aire: 72.70 ± 0.67%; 
Mann-Whitney U test, U = 2578.50, p < 0.0001, Table 1). 


Daily emergence and retreat 

Lizards emerged from their shelters on Aire and began 
their activity about 30 min earlier than on Colom (Table 2). 
In addition, lizards on Aire spent more time on their first 
basking session before beginning to forage. Both humid¬ 
ity and wind were higher on Aire than on Colom during 
the four study periods. Air temperatures at which activ¬ 
ity started and ended were similar for both islets (Table 2). 


Microhabitat selection 

Lizards from Colom selected rocky microhabitats signifi¬ 
cantly more commonly and avoided microhabitats of grass 
and bushes (G = 89.688, 3 d.f, p < 0.0001; Table 1, Fig. 1). 
They also preferred microhabitats in full sun and actively 
avoided fully shaded microhabitats (G = 30.859, p < 0.0001; 
Table 1). Lizards on Aire selected rocky microhabitats and 
avoided microhabitats of grass and bushes (G = 38.067, 
3 d.f, p < 0.0001; see post-hoc comparisons in Table 1, 


full shade) were similar for the three thermoregulation sit¬ 
uations (‘cold’, ‘optimal’, and ‘warm’) on Colom (G = 7.040, 
4 d.f, p = 0.134), as well as on Aire (G = 2.690, 4 d.f, p = 
0.611; Fig. 2). 

On Golom, a univariate test of abiotic variables demon¬ 
strated that lizards actively selected microhabitats with sig¬ 
nificantly higher substrate and air temperatures (Table 3) 
while wind speed was not significant (Table 3). Microhabi¬ 
tats with higher relative humidity were significantly more 
often avoided by ‘optimal’ and ‘cold’ lizards (Table 3). The 
paired logistic regression model that better explains the 
data of ‘cold’ lizards on Golom islet revealed that substrate 
temperature was the main abiotic trait affecting microhabi¬ 
tat choice, whereas the remaining abiotic traits, air temper¬ 
ature, wind speed and humidity, had no significant impacts 
(Table 4). For ‘optimal’ lizards on Golom, both air and sub¬ 
strate temperatures significantly affected the choice of mi¬ 
crohabitat, whereas wind speed and humidity did not have 
significant effects (Table 5). An increase by 10% in air tem¬ 
perature (approximately 2.5°G) would increase the proba¬ 
bility of selection of a microhabitat by 19%, and an increase 
by 10% (approx. 2.8°G) would increase the probability of 
selection by 23.7%. For ‘warm’ lizards, none of the studied 


Colom island 


Aire island 



Rock Soil Grass Bush Rock Soil Grass Bush 


Figure 1. Proportion (%) of each microhabitat type selected by 
Podarcis lilfordi lizards and their availability on both islets. 


microhabitats under each type of insolation situation (full sun, 
filtered sun, or shade) according to three thermal situations: 
(1) cold (T|^ < PTR [preferred temperature range]), (2) optimal 
(T|^ within the PTR), and (3) warm (T^^ > PTR), for both islets. 
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Table 3. Mean ± SE (n) values of four microclimatic variables on Colom, comparing the capture points of lizards (cases) to the avail¬ 
ability of alternatives (controls) in three thermal situations: (1) cold (Tj^ < PTR), (2) optimal (T^ within PTR), and (3) warm (T^ > PTR). 


Variable 

Situation of 
thermoregulation 

Lizards 

Availability 

F 

P 

Substrate temperature 

Cold 

25.55±0.61 (41) 

23.55±0.30 (155) 

9.245 

0.003 


Optimal 

29.25±0.54 (53) 

27.25±0.37 (155) 

8.104 

0.005 


Warm 

31.52±0.50 (37) 

29.12±0.40 (92) 

11.599 

0.001 

Air temperature 

Cold 

23.43±0.47 (41) 

22.13±0.21 (155) 

7.610 

0.006 


Optimal 

26.86±0.78 (53) 

24.63±0.22 (155) 

22.519 

< 0.0001 


Warm 

28.66±0.46 (37) 

25.69±0.30 (92) 

28.218 

< 0.0001 

Wind speed 

Cold 

0.24±0.05 (41) 

0.33±0.03 (155) 

3.351 

0.069 


Optimal 

0.43±0.06 (53) 

0.44±0.23 (155) 

0.185 

0.668 


Warm 

0.56±0.06 (37) 

0.48±0.04 (92) 

0.150 

0.699 

Relative humidity 

Cold 

70.63±1.70 (41) 

70.64±0.77 (155) 

0.000 

0.988 


Optimal 

59.19±1.19 (53) 

61.75±0.63 (155) 

3.994 

0.047 


Warm 

55.03±1.63 (37) 

58.80±0.69 (92) 

6.324 

0.013 


Table 4. Mean ± SE (n) values of four microclimatic variables on Aire, comparing the capture points of lizards (cases) to the availability 
of alternatives (controls) in three thermal situations: (1) cold (Tj^ < PTR), (2) optimal (Tj^ within PTR), and (3) warm (Tj^ > PTR). 


Variable 

Situation of 
thermoregulation 

Lizards 

Availability 

F 

P 

Substrate temperature 

Cold 

26.68±0.68 (22) 

24.78±0.30 (99) 

7.687 

0.006 


Optimal 

30.19±0.49 (73) 

28.09±0.30 (290) 

10.751 

0.001 


Warm 

32.69±0.99 (11) 

30.76±0.61 (44) 

2.154 

0.148 

Air temperature 

Cold 

24.85±0.44 (22) 

23.73±0.19 (99) 

6.977 

0.009 


Optimal 

26.86±0.37 (73) 

25.27±0.16 (290) 

17.469 

< 0.0001 


Warm 

28.58±0.55 (11) 

26.47±0.28 (44) 

11.204 

0.002 

Wind speed 

Cold 

0.58±0.13 (22) 

0.56±0.06 (99) 

0.004 

0.953 


Optimal 

0.62±0.72 (73) 

0.55±0.04 (290) 

0.664 

0.416 


Warm 

0.64±0.21 (11) 

0.59±0.10 (44) 

0.044 

0.834 

Relative humidity 

Cold 

74.23±2.23 (22) 

75.06±0.68 (99) 

0.572 

0.451 


Optimal 

70.38±1.09 (73) 

72.74±0.44 (290) 

5.153 

0.024 


Warm 

66.36±1.91 (11) 

68.89±0.76 (44) 

1.970 

0.166 


abiotic constraints had a significant effect on the micro¬ 
habitat choice offl lilfordi brauni lizards (Table 5). 

On Aire, univariate tests demonstrated that lizards pre¬ 
ferred microhabitats with significantly higher air tempera¬ 
tures, regardless of their body temperatures (Table 4). How¬ 
ever, only ‘cold’ and ‘optimaF lizards selected microhabitats 
with higher substrate temperatures (Table 4). Wind speed, 
again made no difference for lizard microhabitat choice, 
and microhabitats with high humidity were avoided only at 
‘optimal’ body temperatures (Table 4). The paired logistic 
regression model that better explains data of ‘cold’ lizards 
on Aire reveals that none of the studied abiotic constraints 
had a significant effect on lizard microhabitat preferences 
(Table 6). For ‘optimal’ lizards, air and substrate tempera¬ 
tures as well as wind speed significantly affected the lizards’ 
choice of microhabitat (Table 6). An increase by 10% in air 
temperature (approx. 2.5°C) would increase by 47.5% the 
probability of selection of a microhabitat, while an increase 


by 10% in substrate temperature (approx. 2.9‘’C) would in¬ 
crease the probability of selection by 14.6%, and an increase 
of a 10% in wind speed (approx. 0.06 ms ‘) would increase 
the probability of selection of a microhabitat by 79.2%. Fi¬ 
nally, the paired logistic regression model that better fits 
data (lowest AIC^) of ‘warm’ lizards on Aire, demonstrated 
that none of the studied abiotic constraints had a signifi¬ 
cant effect on lizards’ microhabitat choice (Table 6). 


Discussion 

Body temperatures of P. lilfordi during spring of 2007 were 
similar on both islets, with a mean T|^ of 32.85‘’C. These 
results differ from those of a previous study in 2006, in 
which Balearic lizards showed a mean T^^ of approximate¬ 
ly i°C higher on Colom during the same period (Perez- 
Mellado et al. 2013, Ortega et al. 2014). A similar situ- 
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Table 5. Paired logistic regression models that best predict microhabitat selection of Podarcis lilfordi on Colom. The model cold’ analy¬ 
ses microhabitat selection of lizards with body temperatures (T^^) below the preferred temperature range (PTR). The model optimal’ 
analyses microhabitat selection of lizards in which Tj^ fall within the PTR, and the model of ‘warm’ analyses microhabitat selection 
of lizards with Tj^ higher than the PTR. 


Model 

AIC 

c 

Variable 

Coeff. 

SE 

P 

Odds ratio 
(10% change) 

95% Cl odds ratio 
(10% change) 

Cold 

29.3 

relative humidity 

0.28 

0.17 

0.098 

1.32 

(0.950, 1.838) 



substrate temperature 

1.43 

0.50 

0.004 

4.18 

(1.561, 11.202) 



wind speed 

-3.43 

2.17 

0.114 

0.03 

(0.001, 2.291) 

Optimal 

27.4 

relative humidity 

-0.19 

0.12 

0.127 

0.83 

(0.652, 1.055) 



air temperature 

0.64 

0.27 

0.016 

1.90 

(1.125, 3.202) 



substrate temperature 

0.86 

0.37 

0.021 

2.37 

(1.141, 4.906) 



wind speed 

-3.57 

1.92 

0.063 

0.03 

(0.001, 1.210) 

Warm 

6.6 

air temperature 

7.478e+01 

4.798e+03 

0.988 

2.988e-i-32 

O-oo 



substrate temperature 

-1.063e+01 

9.861e+02 

0.991 

2.422e-05 

O-oo 



wind speed 

6.282e+02 

4.000e+04 

0.987 

6.406e-i-272 

O-oo 


Table 6. Paired logistic regression models that better predict microhabitat selection of Podarcis lilfordi on Aire. The model ‘cold’ analyses 
the microhabitat selection of lizards with body temperatures (Tj^) below the preferred temperature range (PTR). The model ‘optimal’ 
analyses microhabitat selection of lizards in which Tj^ fall within the PTR, and the model of ‘warm’ analyses microhabitat selection 
of lizards with Tj^ higher than the PTR. 


Model 

AIC 

c 

Variable 

Coeff. 

SE 

P 

Odds ratio 
(10% change) 

95% Cl odds ratio 
(10% change) 

Cold 

18.6 

air temperature 

1.14 

0.68 

0.096 

3.12 

(0.817, 11.929) 



substrate temperature 

0.76 

0.60 

0.206 

2.13 

(0.659, 6.884) 

Optimal 

44.8 

air temperature 

1.56 

0.42 

< 0.001 

4.75 

(2.066, 10.935) 



substrate temperature 

0.38 

0.16 

0.020 

1.46 

(1.061, 2.004) 



wind speed 

2.07 

0.92 

0.024 

7.92 

(1.307, 47.945) 

Warm 

4.6 

air temperature 

3.414e+01 

4.368e+03 

0.994 

6.708e+14 

O-oo 



wind speed 

3.465e+01 

4.601e+03 

0.994 

1.119e+15 

O-oo 


ation was observed for air and substrate temperatures at 
the places where lizards were captured. Although mean 
temperatures at the closest weather station were similar in 
both years (data available from the Spanish Meteorologi¬ 
cal Agency, ‘Agenda Espanola de Meteorologia, AEMET’; 
http://www.aemet.es), it is possible that the spring of 2007 
was harsher on Colom than on Aire. These results empha¬ 
size the thermal variability to which Balearic lizards are ex¬ 
posed from one islet to another. Wind at the capture places 
was twice as fast on Aire than on Colom, while the relative 
humidity was also significantly higher on Aire. It is possi¬ 
ble that a higher relative humidity may be related to faster 
wind speeds, since the wind blows directly from the sea 
on these little islets. Either way, the wind speed would be 
the main abiotic difference between both islets during the 
study period. Mean abiotic conditions of the habitat, meas¬ 
ured at the four points surrounding each lizard’s capture 
point, both air and substrate temperatures, as well as wind 
speed and relative humidity were higher on Aire than on 
Colom. 

An interesting result is that lizards on Aire emerged 
from their shelters approximately half an hour before liz¬ 


ards on Colom and took longer basking sessions before 
they started activity. They emerged and started their activ¬ 
ity when the air temperature was approximately 17-19°C. 
They also retreated to their shelters more than half an hour 
later on Aire than on Colom, so that lizards on Aire have a 
longer diel activity than lizards on Colom. This result may 
be related to the melanistic coloration of Aire lizards, as 
has been suggested in previous studies (Perez-Mellado 
et al. 2013, Ortega et al. 2014). A tailor-made experimental 
study focusing on cooling and heating rates of lizard from 
both populations would clarify this question. Nonetheless, 
both populations ended their activity periods at approxi¬ 
mately 20°C and retreated to their shelters at approximate¬ 
ly i8°C. Differences in activity as a means of thermoregula¬ 
tion are common, and have been reported for other lacer- 
tids as a way of dealing with habitat variations, for example, 
in altitude (DIaz 1997, GvozdIk 2002). A longer diel activ¬ 
ity period could compensate for the stronger winds on Aire 
and also allow improving lizard physiological performance 
in this population (e.g., Adolph & Porter 1993). Although 
melanistic coloration has been reported to play a role in in¬ 
creasing the fitness of some ectotherms living in cooler ar- 
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eas (Clusella-Trullas et al. 2007), we still do not know 
if the melanistic coloration is related to these thermal ad¬ 
vantages in the Aire population. It is possible that main¬ 
taining high body temperatures for longer periods of activ¬ 
ity would allow continuous active foraging in P. 1 . lilfordi, 
and thus would facilitate the highly plastic omnivory of the 
population (Perez-Mellado 1989, Perez-Mellado & 
CoRTi 1993, Perez-Cembranos & Perez-Mellado 2015, 
Perez-Cembranos et al. 2016) in a habitat that is almost 
free of predators (e.g., Cooper & Perez-Mellado 2004, 
MencIa et al. 2017). In fact, one of the characteristics of 
small isolated islets is that they have limited food resourc¬ 
es. Consequently, foraging activity has to be extended for 
longer periods to satisfy the daily energetic requirements 
of lizards (Brown & Perez-Mellado 1994). 

Balearic lizards selected sunny microhabitats and avoid¬ 
ed shaded places, similar to other lacertids (e.g., DIaz 1997, 
ScHEERS & VAN Damme 2002). Shuttling between sun and 
shade is also a common thermoregulation mechanism 
(DfAZ 1997, Castilla et al. 1999). We report a similar pat¬ 
tern of use of sun/shade patches from the Balearic lizards 
from both islets. This pattern could also be similar to that of 
other Mediterranean lacertids. For example, shuttling be¬ 
tween sunny and shaded patches in the morning is an im¬ 
portant mechanism in P. liolepis on the Columbretes Islands 
to avoid overheating over midday in summer (Bauwens et 
al. 1996). Furthermore, Balearic lizards selected rocky mi¬ 
crohabitats and avoided grass microhabitats, especially on 
Colom, where grassy areas were highly available. The selec¬ 
tion of rocky substrates would facilitate faster heating rates 
(Beillure & Carrascal 2002, Carrascal et al. 2002). 
More research on the behavioural thermoregulation of Bal¬ 
earic lizards would help us to fully understand the mecha¬ 
nisms of thermoregulation on these variable little isles. 

Interpreting microhabitat selection in relation to ther¬ 
mal requirements was confusing on Aire, because only liz¬ 
ards with optimal body temperatures (31.78-35.68°C) ap¬ 
peared to select microhabitats based on their abiotic traits. 
In contrast, lizards on Aire with body temperatures low¬ 
er or higher than the preferred range did not choose their 
microhabitats based on their thermal characteristics. This 
result would suggest that other factors, like food availabil¬ 
ity, could influence the selection of microhabitats on Aire 
during spring. On Colom, P. lilfordi selected warmer mi¬ 
crohabitats from those available when their body tempera¬ 
tures were below the optimal range. Thus, Colom lizards 
selected their microhabitats for warming according to their 
substrate temperatures, while air temperature, wind speed 
and humidity did not affect these choices. This is congru¬ 
ent with previous results from other species, since lizards 
prefer the substrates that facilitate faster warming rates 
(Beillure & Carrascal 2002). In addition, when lizards 
on Colom had optimal body temperatures, they also pre¬ 
ferred the microhabitats that offered warmer substrates 
and higher air temperatures. Finally, lizards on Colom with 
body temperatures above their optimal range did not show 
any preferences for microhabitats according to their ther¬ 
mal traits. 


In conclusion, wind and humidity did not affect micro¬ 
habitat selection of the two Balearic lizards studied, while 
substrate temperature was the main factor driving the mi¬ 
crohabitat choice of these lizards. Although a stronger 
wind may potentially have a cooling effect on substrates 
and, thus, affect thermoregulation, the wind speeds on 
Aire (0.61 ms 0 were two times higher than on Colom 
(0.29 ms ‘) and did not affect the lizards’ microhabitat se¬ 
lection. Logan et al. (2015) found that wind speed affected 
the activity of Anolis lemurinus, whereas humidity did not. 
We also found a detrimental effect of wind on thermoregu¬ 
lation effectiveness in Iberolacerta aurelioi, when the wind 
freshened up from 0.71 to 1.73 ms ‘ and so potentiated the 
cooling of substrates (Ortega et al. 2017). Since climate 
and habitat interact in shaping the thermal reaction norms 
of lizards (Rutschmann et al. 2016), a next step could be 
to study the roles of the abiotic factors that affect micro¬ 
habitat selection in shaping thermal reaction norms in the 
Balearic lizard. 
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